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Abstract 


We consider the eigenvalues of the fractional-order 
Sturm-—Liouville equation of the form 


—©Doy o Dor y(t) + aul) = Ay(t), O<a<l, te [0,1], 
with Dirichlet boundary conditions 
Ihe °y@leso =0 and I}y*y()|e=1 = 0, 
where q € L?(0,1) is a real-valued potential function. The method is used 


based on Picard’s iteration procedure. We show that the eigenvalues are 
obtained from the zeros of the Mittag-Leffler function and its derivatives. 
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1 Introduction 


Fractional Sturm—Liouville Problems (FSLPs) are generalizations of the clas- 
sical Sturm—Liouville Problems in which the ordinary derivatives are replaced 
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by fractional derivatives or derivatives of fractional order. As an introduc- 
tion, the interested reader may wish to consult the great variety of works on 
the subject, starting with books such as [14, 20, 23]. Several authors have 
considered the numerical FSLP, for example Al-Mdallal [2, 24] applied the 
Adomian decomposition method for solving fractional Sturm—Liouville prob- 
lems. Abbasbandy and Shirzadi [1] applied the homotopy analysis method 
for solving fractional Sturm—Liouville problems. Also in [8], the eigenvalue 
problems for the fractional ordinary differential equations have been inves- 
tigated with different classes of boundary conditions including the Dirichlet, 
Neumann, and so on. They explained that choosing a = 2 leads to the 
classical ones of the second-order differential equations. When the order a 
satisfies 1 < a < 2, the eigenvalues can be finitely many; see [9, 22]. It has 
been applied to many fields in science and engineering, such as viscoelas- 
ticity, anomalous diffusion, fluid mechanics, biology, chemistry, acoustics, 
control theory, and so on. In the applications mentioned above, a class of 
integro-differential equations with singularities, fractional differential equa- 
tions have been involved; see [11, 16, 15, 17, 21, 18, 20, 5]. In [25], the 
authors have considered a regular fractional Sturm—Liouville problem (RF- 
SLP) of two kinds RFSLP-I and RFSLP-II of order v € (0,2) with the 
fractional differential operators both of Riemann-—Liouville and Caputo type, 
of the same fractional-order 4 = v/2 € (0,1). It was proved that the reg- 
ular boundary-value problems RFSLP-I & -II are indeed asymptotic cases 
for the singular counterparts SFSLP-I & -II. The inverse Laplace transform 
method for obtaining analytical solutions of the FSLPs and eigenvalues has 
been investigated in [10]. The reproducing kernel method has also been used 
to calculate the eigenvalues of the FSLPs. Dehgan and Mingarelli [6, 7] have 
investigated the general solution of three- or two-term fractional differential 
equations of mixed Caputo/Riemann-—Liouville type in the case of Dirichlet 
boundary conditions. From a numerical viewpoint, we also refer the reader 
for fractional differential equations to [2, 3, 4, 12, 24]. Particularly, the 
boundary value problem is of the form 


— D2. o Da.y(t) + q()y(t) = Aw), O<a<l, t€ (01, (2) 


Le ula =0,; and Tyo =0. 


For each 1/2 < a < 1, it is proved that there is a finite number of real 
eigenvalues, an infinite number of nonreal eigenvalues, that the number of 
such real eigenvalues grow without bound as a + 1~. Also, they expressed 
asymptotically behavior of the eigenvalues as a function of a@ in the following 


form: 
nT 


sin(55) 
This corresponds exactly with the well-known classical asymptotic estimate 
An > n?n? as n — 00; see [6, 7]. The remaining structure of this paper is or- 
ganized as follows: In the next section, we review the Mittag-Leffler function 


oa 


An(a) ~ ( 3 al. 
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and Laplace transform of the Riemann—Liouville and the Caputo fractional 
derivatives. In section 3, we discuss the zeros of the iterative method with 
two examples. In section 4, we discuss the analysis of the iterative method. 
The last section includes our conclusions. 


2 Preliminaries 


We recall some definitions in fractional calculus. We refer the reader to [6] 
for further details. 


Definition 1. The left and the right Riemann—Liouville fractional integrals 
I®, and If. of order a € R* are defined by 


TEL = Bey f -yUFO)as, (a (2) 
and ‘ F 

Te FO) = By f O° Fe)ds, (ad), 
respectively. 


Here ['(a) denotes Euler’s gamma function. The following property is 
easily verified. 


Lemma 1. For a constant C, we have 1%, C = batty -C. 


Definition 2. The left and the right Caputo fractional derivatives “D®, and 
“D,- are defined by 


t 
"Di f(t) = T° o Df) = FW a / (f= a)? -1 7 (aids, i> a, 
(3) 


and for t < b, 


(—1)? 
T(n-—a 


b 
De f(t) = (-L)"o* oD" f(t) = 5 f(s 0 FF (o)ds, (4 


respectively, where f is sufficiently differentiable and n-1l<a<_n. 


Definition 3. The left and the right Riemann—Liouville fractional derivatives 
D®. and D; are defined by 


1 d” 


Dee (t) <= D” oT f(t) _ Tin - a) dt” 


fe s)"-°! F(s)ds, t >a, 
(5) 
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and for t < 8, 


De f(t) = (1D omg (t) = LA" 


b 
Baya | OOO Fl)ds, (6) 


respectively, where f is sufficiently differentiable and n-1l<a<_n. 


2.1 The Mittag-Leffler function 


The function F(z) is defined by 


which was introduced by Mittag-Leffler [20]. In particular, when a = 1 and 
a = 2, we have 

Fi(z) =e*, E2(z) = cosh(V/z). 
The generalized Mittag-Leffler function F4,8(z) is defined by 


k 


Ea,p(2) = > eran (z,0 €C, R(a) > 0), (8) 


of course, when 6 = 1, Eo,e(z) coincides with the Mittag-Leffler function 
(18): 
Eua(z) = Eq(z). 


Two other particular cases of (3) are as follows: 


EF, 2(z) => — = E2(z) = ae. 


Further properties of this special function may be found in [11]. 


Theorem 1 (see [22]). If a < 2, 8 is an arbitrary real number, ju is such 
that 
7 <p <min{z, 7a}, and c is a real constant, then 


|Eo,a(2)| < (u < jarg(z)| <x), lz) 20, z€C. 


1+ |2| 
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2.2 Laplace transform 


Definition 4 (see [22, 19]). The Laplace transform of a function f(t) defined 
for all real numbers t > 0, t stands for the time, is the function F(s) that is 
a unilateral transform defined by 


F(s) = e(7)} = f * em" F(t)at, 


where s is the frequency parameter. 


Definition 5 (see [22, 19]). The convolution of f(t) and g(t) supported on 
only [0, 00) is defined by 


(f*g)(t) = [ f(sg(t—s)ds, fg: 0,00) +R. 


Property 1 (see [22]). The Laplace transform of the convolution of f(t) and 
g(t) is given by following relation: 


L{(f * g)()} = AFM} x Ltgt- 


Property 2 (see [6, 22]). The Laplace transform of the derivatives of the 
Mittag-Leffler function reads as follows: 


- —styak+pB— k a Rist F 1 
| e 8teoktB 1p) (4d )dt = (s® = FET (Re(p) > |a]~). 


Property 3 (see [13, 22]). The Laplace transform of the Riemann—Liouville 
fractional derivative is obtained as 


n-1 
(LoDFy)(s) = s*(Ly)(s) — $5 s”-*- D¥ (ol? *y)(0), n-1<asnneNn. 
k=0 


If0<a<1, then ({oD?y)(s) = s*(Ly)(s) — (off “y)(0). 


Property 4 (see [13, 22]). The Laplace transform of the Caputo fractional 
derivative is obtained as 


(Lp DFy)(s) = s*(Ly)(s) — 3 s°-*-l (pD¥y)(0), n-1l<a<nneENn. 
k=0 


If 0 <a <1, then (£5D¢y)(s) = s*(Ly)(s) — s*~*y(0). 
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3 Iterative method 


If we take the Laplace transformation of equation (1), then we get the fol- 
lowing equation: 


a a-l1 
$ l-a s a 1 
ul) = Spalax *yle=o + Spa Por ul)li=0 + Soa 


L{h(t)}, 


where h(t) := q(t)y(t). We use the inverse Laplace transform 


y(t) =¢,t°~| Baeo(—-AP*) + cat” Fae,a+1(—At?”) 


9 
pO a alae) 40); (9) 


in which, specifying the constants in (2) by setting, without loss of generality, 
c, = 0 and cz = 1 are given constants and * is the convolution symbol. We 
consider the following iterative method: 


Ult) =t" Hin gia (ae) +f (t — 8)?°—" Fou 20(—A(t — 8)?“)q(8)Ym_—1(s)ds. 


Example 1. For simplicity, first we assume that q(s) = 1 and that 
yo(t) = t” Ban, o¢1(—At?”). 


Then we get 


y(t) =t" Fee,a+1(—At?*) 


t 
+ | (t = Bye Bevan —Ab = 8)°*) 8” Fag. o+1(—As?)ds. 
) 
To calculate the integral term, we apply the following method: 


ere f (65)°°—" Boa ae(— AE = 6)" 9" Baga (—As™)ds}) 


= A 430 ppt) 


1! eee Gosia F 


Thus i 
y(t) = t° Boa,ati(—AO*) + Ft Boe a4 (—AP*). 
Now, we calculate y2(t) as follows: 


1 


a p2 a 
9! ° BD (re? ); 


yo(t) = y(t) + 


as the same way, we get 
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(t) = yo(t) + at? BO (X29) 
¥3 = Y2 3) 2a,a+1 ‘ 


Finally, after the iteration, we have the following relation: 


m+1 
1 a ple m 
nosy oo BY sae): (10) 
k=0 ~ 


By choosing some terms from the approximate solution and calculating their 
zeros, the same results are obtained from [7] as follows: 


1 
ya(t) =t" Baaati(—N™) + Ft BS ay (—At™) 


1 a pe a 1 a (3 a 
eal Be ne ea | ee aed 


For example, if we take three terms with the boundary conditions 


Re*uOher = EEN | 9) Puls)as = 0. 


then we get 


_ 1 : wd _ 
euler 2 a may ff tS" Boo wads Vas 


1 ' il a ie 
| rasa | (t — 8) 58 Bae ati (—As"*)ds 


1 t 
+ 7 (t—a)-*5 + 8482) (— \s2%)ds 


1 : il cy, . 
| 5 f= 9) 59 Baa (As) ds} = 0. 


Using the inverse Laplace transform implies 


g2a—-2 


-1 
aah oe pay! 


ies “y3(t)|r=1 = a4 
g2a-2 = 2a—2 
(s2@ Fayit £ {Ta papait 


= [BB 2a.a(—7) + 1041 BE) 9(— Ae) + tt BL 9(— a0) 


+ £-T{ = 


+ OB a(-AP)] = 0. 


Therefore, the zeros of the above relation yield eigenvalues. 


For the case of q(s) = 0, we refer the reader to [6, 7]. 
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Table 1: The eigenvalues \n of the FSLP of Example 1 


koa 0.88 0.92 0.96 0.98 0.99 1 
10. An 10.71568699 10.55992405 10.66410923 10.78678886 10.86333576  10.94942560 
27.42039814  31.48842150 35.71275308 38.01570185 39.23145940 40.49542871 
39.66187430 51.88983922 67.66694524 77.14409122 82.33647923 —87.85647119 
20) An 10.71568686 10.55992405 10.66410923 10.78678886 10.86333576 10.94942560 
27.51414362  31.51192899 = 35.71778637 38.01801903 39.23303636 40.49650499 
60.83896648  66.90351096 76.82733019  82.94041507 86.28709241  89.83417215 
86.83800089  108.6657324 131.4598304 144.4575545 151.4897515 = 158.9183593 
107.4785323  147.2884927 199.9559706 231.8456925 238.3472917 249.1565859 
260.3905771 286.9178097 
275.8685836  300.6931093 


Figure 1: The curves of eigenfunctions, k = 10 for n = 1 (solid line), nm = 2 (dash 
dot line), n = 3 (dash line), where a = 0.98, A1 = 10.78678886, Az = 38.01570185, and 
A3 = 77.14409122 for Example 1. 


Figure 2: The curves of eigenfunctions, k = 10 for n = 1 (solid line), n = 2 (dash 
dot line), n = 3 (dash line), where a = 0.99, Ay = 10.86333576, Ag = 39.23145940, and 
A3 = 82.33647923 for Example 1. 
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Remark 1. Note that as a approaches 1, one can see that the eigenvalues 
satisfy A, = n?n? +1. This shows that our results are a generalization of the 
classical ones. 


Example 2. Now, we assume q(s) = s°, and we have the following iterative 
method, analogous to the previous computation: 


Ym(t) = t* Eoo,0-41(—At?®) +f (t— 8)?" Boq,2a(—X(t — s)?*)q(8)Ym—1(s)ds. 


We define 
yi (t) =t% Boa,a+1(—At?*) 
t 
+f (t oS 8)7°—" Bog 20(—A(t = s)-*) : s°t? Fo at1(—As?”) ds, 
0 
and again to calculate the integral term, we apply the following method: 


£71 (2f 9 Bag,2a(—)} LF? Baa at (Al) }) 


T((2ky +1) + 8B +1) 
— x ky 3 7 (2kh1+3)a+8 pe 
ON eer crs re 2a 


(2k: +3)0-46-41(—At?®). 


CO 


| , T(Qk+la+6+1) 
pea T((2k, + a +1) 


x p2mt3e+s pr, 


(2k +3)a+6+1(—At*%), 


now, we calculate y2(t) as follows: 


yo(t) =yr(t)+ S> So (-AyE te 


k,1=0 k2=0 
T'((2k; +1lha+64+1) T((2ko + 2k, +3)a+ 28 +1) 
T(Qk+la+1)  T((2ke+2k,+3)a+B+4+1) 

x pPhe tobi tote pe 


\(Qha-+2ky+5)0-+26-+1(—At?*), 


(11) 


similarly, we have 
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y3(t) =ye(t) 
ea . I((2k+la+ 8 +1) 
=X kitketks , 
= D4 Dy) T'((Qkh + Na+) 


ky=0 kg=0 k3=0 
T'((2k2 t 2ky t 3)a 26 t 1) ” T((2k3 t 2ke 2k 5)a 30 1) 
T'((2k2 + 2k1 + 3)a+64+1)  T((2k3 + 2ke + 2k + 5)a 4+ 26 4+ 1) 


2k3+2k2+2k1+7)a+3, 2 
pp Net AB OSE Bs (as done vans rhepapta(—At °). 


Finally, we get the following relation for y(t): 


Ym(t) =Ym—1(t) + > oe oy: ki tkot+--+km 


k,=0 k2=0 


ye op Bhim + 2km 1 +2ki4 an 1))a-+mB (12) 


X Ee, (Qk + 2k att 2ky+(2m—1))a+mp-+1)(—At?%), 


where 


A =T((2k1 + l)a + B + 1)T((2ke + 2ky + 3)a + 26 + 1) x 
DT ((2kimn + 2km—1 +++: + 2k, + (2m —1))a+mB +1), 


and 


B=I((2k1 + la + 1)0((2ko + 2k + 3)a+6B4+1)x--- 
x D((2km + 2km-1 +++: + 2k, + (2m —1))a+(m—1)8 +1). 


Now, in order to obtain eigenvalues, by choosing three terms from (37) and 
imposing the following boundary condition, we have 


Thx" ys(t)le=1 


a in (2 aa (t — s)~% 8 Fea,a41(—As?™)ds 
a T'((2k1 +1)a+8+1) 
> ial T'((2k, + 1)a + 1) 


k1=0 


2 ‘ = 2ky+3)at+ 2 
‘. rama (E— 8) 8PM FDO+F Bo, (ony ta)arpti(—As"*)ds 


+ 3 y( yyfathe P(2ki + Dot B41) | D((2ke + 2k + 3)a + 26 +1) 


i: T'((2k1 +1l)a + 1) T'((2ke + 2ky + 3)a+ 641) 
1=0 ko=0 


1 i _ 4 
x a! (t—s) ghar Bhat) e 28 Boa, (2ho+2k1 +5)0-+28-41(—-AS?™)ds 
T(1—a) Jo 
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ai _T((2k1 +1)a+ 8 +1) 
+ —))kitkatks 
, & 2. P((2k1 + La + 1) 


T((2ke + 2k1 +3)a+28+4+1) P((2ks 


TP ((2ko +2k1 +3)a+B+1) T((2k3 


t 
- 1 / (t — 5) 7% s(2hst2k2t2hi +7)a+38 
T(l-a) Jo 


2ko 
2ko 


2ky 
2ky 


5)a+ 3641) 
5)a +2641) 


2 
X Foq,(2k3+2ko+2k1 +7)0+36-41(—AS~°)ds 


3) is os 


Now, by using the inverse Laplace transform, we have 


Ty “ys(t)le=a 


_ 20, \ ne, (2k +)e+8+1) 
tE2a,2(—At?*) 4 2 d) Ge Gea 


x GOR TAOTET By con, +0)0+0+0(—At™) 


a keythy E((2k1 + Ia+B+1) D((2ka + 2k + 3)a + 28 +1) 
7 py Ds a T((2k1 + 1)a +1) T((2k2 + 2ki +3)a+6+1) 


ki =0 kg=0 
x phot2ki t4)ja+26+1 pp 


(Qk +2ky+4)a+28-42(—At?*) 


i: 5 S 3 (—a)hitkeths T((2k1 + l)a+ 8 +1) 
fee cee reah P((2k + la +1) 
T'((2k2 + 2ki +3)a+28+4+1) T((2k3 + 2k + 2k1 +5)a 
T((Qk2 + 2k1 +3)a+8+1) [((2k3 + 2k2 + 2k, +5)a 
7 p2hs+2ko+2k1 +6)a+36+1 pp 


38 
28 


1 
T 
1 
T 


Therefore, solving the above relation with respect to A yields eigenvalues. 


Table 2: The eigenvalues \n of the FSLP for Example 2 


ki a 0.88 0.92 0.96 0.98 0.99 1.0 

ki,ka,k3 =5 = 2An —10.05385926 9.933778015 10.06348862 10.19660718  10.27791426 10.36849015 
47.56007387  32.34315896 35.63420589 37.76603821 38.92583141 40.14691848 

1126.335623 44.07433976 60.04931462 68.65397481 73.15018892 77.78582381 

2383.852442 60.98311732 79.07639037 90.47082662 96.86784919 103.7750278 

4084.543198 1672.025020 2492.533591 3047.891192 3371.626995 3730.657415 

29.99913622 3669.500947 5672.277915 7063.067903 7884.473754 8803.558665 

6468.168446 10283.17818 12984.46721 14595.72336 16410.73873 

44.95131154 62.05287730 72.95773310 79.12198365 85.81612846 

ky, k2,k3 =10 An  10.05563487 —9.934078332 10.06355202 10.19663792 10.27793590 —10.36850547 
26.86376328 30.94976951 —35.18620831 37.49382726 38.71166129 39.97755111 

38.80996910 51.04646579 66.84154959 76.33005616 81.52855664 87.05489991 
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0 0.2 0.4 0.6 0.8 1 


Figure 3: The curves of eigenfunctions, k1,k2,k3 = 5 for n = 1 (solid line), n = 2 (dash 
dot line), n = 3 (dash line), where a = 0.98, A1 = 10.19660718, Az = 68.65397481, and 


A3 = 90.47082662 for Example 2. 
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t 


Figure 4: The curves of eigenfunctions, k1,k2,k3 = 5 for n = 1 (solid line), n = 2 (dash 
dot line), n = 3 (dash line), where a = 0.99, A1 = 10.27791426, Az = 73.15018892, and 


A3 = 96.86784919 for Example 2. 


4 Analysis of the iterative method 


Lemma 2. The n-fold series in relation (37) is convergent fora = 1, 6 EN, 


and |t| < . 
|A| 


Proof. It is sufficient to show that the double series in relation (11) is con- 
vergent. By Theorem 1, we have 


c love) k ! 
lyo(t) — (8) < opty Dao I! SS 


- Qk + 2k +3426)! op (14) 
Ne ( p2kot2ki +5428 
xd il ([(2ky + 2k; +3)) + B)! 


ko=0 
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= |A|, by the basic root test of 


Since lim sup 
kg-00 


k (2k + 2k, +34 28)!|A|F2 
([(2k2 + 2h1 + 3)] + 8)! 
convergence of nonnegative series, we conclude that for |t| < Jor the series 
3 [AIF (2k2 + 2ki + 3 + 28)! p2kot+2k1 +5426 
([2k2 + 2k; +3] + 6)! 


ko=0 
is absolutely convergent. Moreover, it is uniformly convergent on any com- 
act subset of the interval (-—, —-_). 
P Sera, 
On the other hand, the series )7;>_, |A|*2.¢2h2 +241 +5428 is absolutely con- 
vergent on the interval (~~, +) and is uniformly convergent on an 
8 ( Vial Tp y 8 y 
compact subsct of the interval (-—+_, -_). So we can take derivative of it 
e (Jor vin? 
on this interval, term by term, for 6 times and get 


co 


S- ([2k2 + 2k1 + 3H 6 + 1)([2k2 + 2ky +3] + B +2) 
k2=0 


.. ((2k2 + 2ky + 3] + 28)|A|R2e2h2t 21 +5428 


= (8) 
— al a Oa ag, 
ka=0 


= (8) 
= 8 (ancerae > ([Ale?)*) 
k2=0 
p2k1+5+28 


7 o(e) 


Now, using (14), we obtain 


|yo(t ~ml)l 
c hea Sa Age Ve 
Tepe > I +1 ce De ) 


3 ym ee p2ki+5+28 (8) 
1+ |A\t? fet (2k; + 1]! 1 — |A|t? 


m 

B 

r 

ct? p28+4 gp, (2k1 +14 6)! on, 44 a 
= ie :( = > Oi a ) 

B 

A\t? 


< 
a 


ct 


i+ 1—[\ [2k + 1]! 


T+ 


428+4 ™ _— (8) 
(ea Se (2k, +2)--- (2k, +14+ A)t m7) 
ki=0 
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oo B 
- ct? t28+4 yi |A|e2 Garey”) wm 
1+ Alt? (1 — Ale? fm 
= 
oo B 
_ ct? p2b+4 jymegrerer0y® (B) 
14+ |Alt? \1— Alt? 
k,=0 
__ a? 428+4 Cabs rman) - 
14 |Alt? \1— Alt? = 
1= 
ct? p28+4 p6+1 (8) (8) 
~ 1+ vate ye Ge) ) 


Therefore, (14) is appeared to be the derivative of a function and hence is 
convergent. 


2+2k1 +3428)! |A]F2 


It follows from the root test that z = lim sup, +00 4) oA 
|\|. Hence 


(2k2+2k,+3)]+£)! 


C= S5 ([(2kz + 2k + 3)] + 8 + 1)([(2kz + 2k: + 3)] + 8 +2) 


ko=0 
Kee X ([(2k2 oe Qk ae 3)] +4 DB Alh2 ghateee ere 


= (> (jaye etir4erts28) © 


k2=0 


where k = (2k2 + 2k; + 5) — [(2ko + 2k; +5) + 6B — 2]. Then 


= (k) f2k1+5+26 \ (k) 
‘eae (aa \t2 2) = ( ) 
Yale) ae 


k2=0 


2ky +542 
t aii 
A|t? 


Noting that (4 *) is convergent for a= 1 and 6 EN, then 


" 2ki +1 | ¢2k1+5+28\ (k) 
<—o(dprs ide! 
1+ |Alt? \ A 0 [2k +1) + 1]! 1 — Ale? 
Cc pe+2 oo 
S eee 
~ 14 [Ae 6 = |rlt 2 ((e2ks +1)+1])+1)x 


x 


(k) 
(2k +1+ je 


(k) 
Cc to+2 i 
< fe+3 t2k1 (1) 
<orine(one ds ) , 


k,=0 


where | = (2k, + 8) — [(2k1 +3) +1]. We have 
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c tet? 1 
“ po+3 (I)y(k) 
<Tepet pe Tw) <% 


the proof is completed. 


The reader can prove the convergence of triple series and so on in the 
same manner. 


5 Conclusion 


In this paper, we have considered the analytical and numerical solutions of 
the FSLPs and eigenvalue problems for the fractional differential equations 
with Dirichlet boundary conditions by using an iterative method. Moreover, 
we proved that the resulting series is convergent. 
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